Effective Vortex Mass from Microscopic Theory 



O 
O 



o 

u 

Oh' 

C/2 



Jung Hoon Han,^'^ June Seo Kim,^ Min Jae Kim,^ and Ping Ao"^ 

'Department of Physics, Sung Kyun Kwan University, Suwon 440-746, Korea 
^CSCMR, Seoul National University, Seoul 151-747, Korea 
'^Department of Mechanical Engineering, University of Washington, Seattle, WA 98195, USA 

(Dated: February 2, 2008) 

We calculate the eflFective mass of a single quantized vortex in the BCS superconductor at finite 
temperature. Based on effective action approach, we arrive at the effective mass of a vortex as 
integral of the spectral function J{ijj) divided by over frequency. The spectral function is given 
in terms of the quantum-mechanical transition elements of the gradient of the Hamiltonian between 
two Bogoliubov-deGennes (BdG) eigenstates. Based on self-consistent numerical diagonalization of 
the BdG equation we find that the effective mass per unit length of vortex at zero temperature is of 
order m(fe/^o)^ (fc/=Fermi momentum, ^o=coherence length), essentially equaling the electron mass 
displaced within the coherence length from the vortex core. Transitions between the core states are 
responsible for most of the mass. The mass reaches a maximum value at T ~ 0.5Tc and decreases 
continuously to zero at Tc. 
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Phase coherence is the defining characteristic of the 
superfluid matter, and vortices - the quantized twist of 
the underlying phase texture - are the unique elementary 
excitations of the condensate. Their ubiquity as well as 
their role in transport and phase transition makes the 
dynamics of an isolated vortex or their array an integral 
part of our understanding of superfluidity. 

The motion of a single vortex may be phrased as a 
Newtonian equation, My{(Pvy/dt^) = F^, where is 
the force acting on the vortex, and My is its effective 
mass. The position of the vortex is represented by rv,. 
The mass of a vortex is entirely effective in nature 
it requires a force to drag a /i/2e flux quantum across a 
type-II superconductor by an external magnet (assuming 
zero mass for the magnet itself) because, in so moving, 
the vortex interacts with the surrounding electrons. The 
virtual transitions among the quasiparticle states caused 
by the vortex motion lead to renormalization of mass, 
which in this case will account for the entire vortex mass. 
Real transitions between the states, on the other hand, 
lead to dissipation. 

The Newtonian description of the vortex motion is 
a limiting case of the more general quantum mechani- 
cal formulation. Quantum-mechanical law of motion is 
most easily derived in the effective action approach In 
essence, one divides the dynamics into that of quasiparti- 
cles and the vortex, then by integrating out the quasipar- 
ticle degrees of freedom, one is left with the effective dy- 
namics of the vortex alone. This way of studying vortex 
dynamics was first conceived in the work of SimanekQ 
and later extensively developed by Ao and Zhu0,|^. In 
this paper, we derive the mass formula from the effective 
action and calculate it using BCS theory for tempera- 
tures < T < Tc. 

The effective action Seff of a single vortex of a BCS 
superconductor, centered at r„(T) at imaginary time t, 
is given byli HQ 
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The spectral function is the quantity 

J(C^) = ^ J (oH Ea-Eb DlfiEahfiE,) 



(1) 



(2) 



Fermi distribution function with energy x is denoted 
f{x). There is an additional action in S(.jj pertain- 
ing to the transverse motion of the vortex, which we 
do not consider herej^. The states a, 6 are the eigen- 
states of the Bogoliubov-deGennes (BdG) Hamiltonian 
Hq in the presence of vortex at r^, with energies Ea 
and E\, respectively. When we assume that the dynam- 
ics is sufficiently local in timeQ, we may approximate 
r^(r)— r.u(T')w(T — r')r.u(T) and write 



The effective vortex mass is the quantity in parenthesis, 
given by 
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f{Ea) - f{E,) 



{Ea~E,f 
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This formula allows an explicit calculation of the vor- 
tex mass using the BCS theory of superconductivity at 
arbitrary temperature. 

Conclusions in the past vary regarding the vortex ef- 
fective mass. One group of theories predicts a "small" 
mass of roughly one electron mass per atomic length 0, 
whereas another group predicts a "large" mass, of order 
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™(fc/Co)^ {kf— Fermi momentum, = coherence length 
at T = 0) 8]. In either case the theories are hmited to 
zero temperature, and finite-temperature generahzation 
does not seem to be straightforward. 

The eigenstates of the BdG equation in the cylindrical 
coordinates is written in the general form(9j 

Mr) = -i=e--e-« f ) (5) 

in the presence of a gap function A(r) = A(r)e^'^ for 
a rectilinear vortex of length L centered at r = 0. The 
radial functions {ua{r),Va{r)) are obtained as eigenfunc- 
tions of the coupled differential equation (r dependence 
in u,v and A is implicit) 



r'^u" + ru'+[{k'^. + 2mE)r'^-\^fi--j ]u = 2TOr^A?; 

r'^v" + rv' + [{k^-2mE)r^-(^fi+^^ ]v = -2mr^Ai£f) 
I 



where fc^ = kj — k"^. Self-consistency requires that 
the r-dependent gap function obey the relation A(r) = 
VJ2a '^a{r)vair)[l — 2f{Ea)] for some choice of the pair- 
ing interaction strength V. The eigenstates are labeled 
by kz, /i, and a, which labels a set of energy levels for a 
given choice of {kz,fi). Uniqueness of the wavefunction 
requires that /i be half-odd integers. Using these eigen- 
states, the transition amplitude needed in evaluating Eq. 
(@J) can be written out in the form —{a\dHo/drv\f]) — 
{xTiy/2)x 




Upper and lower signs correspond to = A*/3 ^ 1- The 
two eigenstates differ by one unit of angular momentum 
due to the gradient operator which connects them. By 
integrating by parts, Eq. Q is equal to (assuming /Zq = 

/i/3 + 1) 
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We assume a hard-wall boundary condition at r = i?. 
One must be careful in treating the boundary term at 
r = i? on the r.h.s. of the above equation, which is 
generally non-zero. 

We follow earlier worksj^ and write the eigenfunc- 
tions in the form 

Ua{r) = ^CQj0„i(r), Va{r) = ^ c?Qi0„+i,i(r), 

i i 

V2 ( r \ 

<l>m {r) = — r Jn I Zm ^ ■ (9) 

RJn+l(Zni) V R/ 

Zni is the i-ih zero of the Bessel function J„. Integer value 



n is related to the angular momentum hy — n -\- i^. 
Coefficients {cai,dai) are determined from matrix diag- 
onalization. Negative-/^ states need not be considered 
separately, instead one can use a positive-^ eigenstate 
{ua,Va) of energy Ea to construct a negative- /i eigen- 
state, given by (wq, — Mq), of opposite energy —Ea- 

With the eigenfunctions thus obtained, we evaluate the 
vortex mass using Eq. and the transition ampli- 

tude given by Eq. (jHJ. Furthermore, Eq. can be 
re-expressed using the coefficients {cai,dai) and making 
use of Bessel identities. For fia = ?i + 5, M/J = ^ ~ 5j 
n > 0, Eq. ||SJ) becomes 
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and for fia — /i^ = — i, it equals 



En — En 
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(10) 



-t>Q/3 — 2^ CaiClpjZniZij 



Mzij) 



2 2 



i?2 



2 , Jo{zii) I Ea + Ep ^0 

UaiCpjZiiZoj 



' J2{zu) \ zf 



Oj 



''il 



i?2 



(11) 



In deriving Eqs. (|10|I -H11 |I we expressed energy in units 



of the zero-temperature gap value at large distance from 
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the core, denoted A(0), and the length in units of kj^. 
In Eq. H10|l . a and jS refer to eigenstates with a given 
angular momenta fia and fip, each. We made use of the 
mapping between negative-/i and positive-/! eigenstates 
to express the result in Eq. solely in terms of the 

fi = +^ eigenstates. The vortex mass can be calculated 
using My = 
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(12) 



Factor 2 multiplying the first term reflects the negative-/* 
contributions which equals those of the positive-/^ tran- 
sitions. 

Equations I|1U|1 - I|12|) complete the formal derivation of 
the vortex mass formula in the superconductor. Using 
Eqs. (|10|I - H11II gives dramatic improvement in both the 
time and accuracy of the calculation of the transition 
amplitude over the brute-force numerical integration of 
Eq. or ©. 

In the following we discuss the vortex mass calculated 
using Eq. H12|l . with the eigenstates obtained from self- 
consistent numerical diagonalization of Eq. ©. We 
choose the coherence length Ef /A{0) = ~ 20 and 
the radius of the boundary R = 100. For each angular 
momentum /i, eigenstates with energies within ±10A(0) 
were retained. This left us with about 90 eigenstates for 
/i = 1/2, and a decreasing number of states for larger 
/i. Calculations were restricted to fcz = only. Vary- 
ing kz will lead to different effective Fermi momentum 
kr = \Jk'j — k1 in Eq. lO and will not change the quali- 
tative conclusion from the k^ — Q case. The energy gap 
vanishes completely at T « 0.571A(0), which we take 
as the transition temperature Tc- We work at several 
(reduced) temperatures t = T/Tc in the range < i < 1. 

Figure ^ shows the calculated gap at large distance 
from the core, A(r), and the /i = +^ core energy level, 
ei(T), for the temperatures we considered. Both quan- 
tities decrease monotonically with T. For T sufficiently 
close to Tc we can no longer resolve the core levels as 
distinct from the continuum. 

Referring to Eq. (|12|l we may divide the mass as aris- 
ing from transitions between the core states [M^'^) , core- 
to-extended states (M^*^), and finally between extended 
states (M„^^), as M„ = M^" + M^" + M^". As it turns 
out the matrix element, Eq. |(7J), is vanishingly small 
between a core and an extended state, and the mass is 
effectively M„ w M^'^ + In Fig. I^^a) we show the 

total mass My(T) and the core-to-core mass, Af^'^(r). 
As is evident from the figure, M^'^iT) > M^''{T) for all 
temperatures except very near Tc where core levels are 
not resolved, but here the total mass is vanishingly small 
anyway. Both the total mass My and the level-resolved 




FIG. 1: Pair potential well away from the core, A(r) (blue), 
and the core energy level for ^ = -|-^, ei/2(r) (red), plotted 
vs. T /Tc- Both quantities are normalized by their respective 
T = values, with ei/2(0)/A(0) = 0.088 from self-consistent 

calculation. A(r) shows excellent fit to A(0)i/l - (T/Tc)^-^ 
(black curve). 



masses M^"^ and M^*^ reach a maximum value around 
T = 0.5Tc. 

The core-level contribution to mass can be further 
grouped according to the angular momentum channels 



n + 



{fjLp = n — i) over which the tran- 



sition takes place. We denote such angular-momentum- 
resolved, core-to-core mass, My'^'"{T). Figure|2Ib) shows 
M^'^"{T) and En^o ^v"'"" (T) = M^%T)- M^''^''{T). We 
find that M^'^'°(r), between = ±i core levels, sur- 
vives at zero temperature and monotonically decreases 
at higher T. The n ^ channels give zero mass for 
T — 0, reaches a maximum at some intermediate T, and 
decreases to zero at T = Tc- This behavior accounts for 
the observed maximum in the mass My(T). 

The results shown in Fig. EJb) can be nicely under- 
stood thanks to the approximate identity that holds for 
the core levels which states | (a| V„i/"o|/3) | ~ \Ea — Ep\. 
We have verified that this relation holds with excellent 
accuracy for all temperatures, except when T ^ Tc. 
Then, we can approximate (see Eq. 1)12(1 ^ 



M-'"(r) « 2meo 



/(e„_i(r))-/(e„+i(T)) 



(13) 



In particular, for n = we have the mass M^'^'°(T) = 
2m^otanh(/3ei (r))/2ei(T). At zero temperature we get 
M='='0(0) = 2mCo/2ei(0) ~ since A(0)/ei(0) ~ 
This is precisely the mass of the electrons occupying the 
area ~ ^q. Using the ei{T) values shown in Fig. we 
find that the approximate formula in Eq. H13|l for n — 
reproduces M^'^'°(T) obtained from the more general 
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FIG. 2: (a) Effective mass M„(T) and the core-to-core con- 
tribution to mass My'^{T) at temperature T/Tc, normalized 
by M„(0). One finds M„(T) « M='=(T) over most of T. (b) 
M^^iT) = MT'^^iT) + E„_^oM^'='"(T) (see text for defini- 
tion). M=='°(r) and M^^iT) - M^^'^iT) are plotted sepa- 
rately. MS'^'°{T) is monotonically decreasing while the higher- 
n channels give a maximum mass at T/Tc ~ 0.5. The core 
levels are not well resolved for T too close to Tc, which ex- 
plains the absence of data points for T ~Tc. 



for a discussion of vortex friction). 

Suppose now tliat tlie core levels got smeared due 
to the impurities. If all higher angular momentum 
core levels except n — Q disappeared, we would expect 
M^'='"5^°(r) is effectively zero, and the mass behaves as 
Mv{T) Ri A/^'^'°(r) showing the monotonic behavior of 
Fig. |2Ib). If all core levels disappeared, wc will be left 
with the extended states contributions which are shown 
to be quite small. It is conceivable that in this case the 
mass becomes "small" , in agreement with predictions of 
Ref. 0. 

In closing we mention that studying the size (R) de- 
pendence may be necessary to get a better estimate of the 
extended-states mass, M^'^{R). Since the core states are 
localized within a few coherence lengths from the origin, 
our choice R — lO^o is presumably adequate to conclude 
M^^iR = 100) « M^'={R = oo). Extended states do not 
meet this requirement, at least for ultra-clean systems 
with no impurity scattering to provide the cutoff. The 
question of friction^, 0, which involves transitions 
between same-energy extended states, will also require 
the study of i?-dependence of the transition rates. We 
will return to these issues in a future publication. 
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formula, Eq. (|12|l . with very high accuracy. 

On the other hand, M^^'^^^iT) w 
/3/(4cosh2(/?e„(r)/2)), e„(r)^^ (e,„+ 1 (T) + e„_ i (r))/2, 
provided the temperature is much larger than the 
typical core energy spacing, T > |e,„+i (T) - e,„_ i (r)|. 
As the typical core energy spacings are a few percent 
of the energy gap A(T), this is not a very restrictive 
condition except near T — 0. Then one can easily 
check that this approximate form for M^'^'"(T) rises to 
a maximum value for T « 0.69e„ provided e„ is treated 
as temperature- independent. Therefore, the maximum 
in the mass shown in Fig. |2Ia) is due to the core-to-core 
transitions in the higher angular momentum channels. 

Summarizing our findings, (i) At zero temperature the 
effective mass My{T = 0) is the mass of the electrons 
forming a cylinder of radius ^o- K increases upon higher 
temperature, reaching a maximum at T « O.STc, and 
vanishes at Tc. (ii) The transition between the localized 
eigenstates forming the core spectrum is mainly responsi- 
ble for the effective mass. Hence, My{T) « M^^iT). On 
the other hand, the dissipation experienced by the mov- 
ing vortex is due to the extended states as it requires 
transition between states of the same energy (See Refs. 
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